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Abstract. Active prominences exhibit plasma motions, resulting in difficul- 
ties with the interpretation of spectroscopic observations. These solar features 
being strongly influenced by the radiation coming from the solar disk, Doppler 
dimming or brightening effects may arise, depending on which lines are observed 
and on the velocity of the plasma. Interlocking between the different atomic en- 
ergy levels and non local thermodynamic equilibrium lead to non-trivial spectral 
line profiles, and this calls for complex numerical modeling of the radiative trans- 
fer in order to understand the observations. We present such a tool, which solves 
the radiative transfer and statistical equilibrium for H, He I, Hell, and CaH in 
moving prominences where radial plasma motions are taking place. It is found 
that for isothermal, isobaric prominence models, the Hell resonance lines are 
very sensitive to the Doppler effect and thus show a strong Doppler dimming. 
The Call lines Doppler effect for the prominence models considered here. We 
illustrate how the code makes it possible to retrieve the plasma thermodynamic 
parameters by comparing computed and observed line profiles of hydrogen and 
helium resonance lines in a quiescent prominence. This new non-LTE radiative 
transfer code including velocities allows us to better understand the formation 
of several lines of importance in prominences, and in conjunction with observa- 
tions, infer the prominence plasma thermodynamic properties and full velocity 
vector. 



1. Introduction 



The main motivation of our modeling work is to contribute to building realistic 
prominence models. For this, we need an accurate knowledge of thermody- 
namic quantities (temperature, densities, ...), level populations (useful, e.g., 
to infer the magnetic field properties from spectro-polarimetric observations), 
velocity fields, energy budget. However these quantities still have large un- 
certainties associated with them. Observations of several different lines from 
different atoms/ions allow us in theory to measure these parameters. Among 
these lines, the H and He lines are important as they are strong and largely 
contribute to the radiative losses. However the prominence plasma being out 
of LTE and optically thick in H and He resonance lines, the interpretation of 
line spectra or intensities in radially moving prominences is a non-trivial task. 
Therefore, non-LTE radiative transfer calculations including velocity fields are 
needed to build realistic prominence models. Here we present such calculations 
and preliminary results. 
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2. Modeling Procedure 

The prominence is represented by a ID plane-parallel slab standing vertically 
above the solar surface. Each prominence model is defined by a set of free pa- 
rameters: the temperature T, the gas pressure p, the slab thickness L (or the 
total column mass), the height of the slab above the limb H, the microtur- 
bulent velocity, and the radial speed. For this preliminary study we consider 
isothermal and isobaric prominences, although the code allows for inclusion of 
a transition region between the cold prominence and the hot corona. We first 
solve the pressure equilibrium, the ionization equilibrium, and the coupled sta- 
tistical equilibrium (SE) and radiative transfer (RT) equations for a 20 levels 
H atom. Then the SE and RT equations are solved for other elements: He I 
(29 levels) and Hell (4 levels), and CaH (5 levels). More details on the mod- 
eling of the hy drogen, calcium, and hel i um spectra in quiescent prorn i nence s 
can be found in Gouttebroze &: Labrossd (|2000l ): lGouttebroze fc Heinzell (l2002l ): 



Labrosse &: Gouttebroze (|2004l ) respectivelv. and references therein. 

For the modeling of active and eruptive prominences, we use a velocity- 
dependent incident radiat ion as boundary condition s for the RT equations. It 



has already been shown bv lHeinzel &: Rompoltl ( 19871 ) in the case of the hydrogen 



lines that the Doppler effect induces a frequency shift of the incident profile 
relative to the rest case, and also a distortion of the incident profile due to 
the variation of the Doppler shift with the direction of the incident radiation. 
It is also the case for the helium ( Labrosse et al. 20061 ) and calcium incident 



radiation. 



3. Diagnostics of Radial Velocities 
3.1. Integrated intensities 

We reproduce the results of iGontikakis et all ( 1997al lbl) who computed the hy- 



drogen radiation emitted by a radially moving prominence, using partial redistri- 
bution in frequency (PRD) for the Lyman lines. We obtain the same variation 
of the relative intensities (intensities normalised to the line intensities when 
the prominence is at rest) and the same line profiles for Lya, Ly/?, and Ha. 
The main result is that there exists an important coupling between Ly/J and 
Ha which causes these lines to be first Doppler brightened, and then Doppler 
dimmed, with increasing velocity, while there is only a Doppler dimming effect 
on Lya. 

Figure [1] presents relative intensities as a function of velocity for the He I 
584 A, Hell 304 A, and He I 10830 A lines (left panel) and Call K and Call 
8542 A lines (right panel) at two different temperatures (8000 K and 15000 K). 

The He I 10830 A line does not show any sensitivity to the Doppler effect, 
which is mainly due to the very weak incident absorption line. The He I 584 A 
line is quite sensitive to the Doppler effect. Its Doppler dimming is more impor- 
tant at low temperature. The He II resonance lines are the most sensitive to the 
radial velocity of the plasma (the relative intensity of the Hell 256 A line, not 
shown, exhibits a similar variation as Hell 304 A), and the Doppler dimming is 
strong at the temperatures considered in this study. Such a result was expected 
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Figure 1. Relative intensities (intensities normalized to the line intensity 
at rest) as a function of velocity at 8000 K (solid lines) and 15000 K (dotted 
lines) for left: He I 10830 A (diamonds), He I 584 A (plus), and Hell 304 A 
(stars); right: Call K (plus) and Call 8542 A (stars). 



since the main mechanism of formation at these temperatures for thes e lines is 
the scattering of the incident radiation ( Labrosse Gouttebroze 200 ll ). Let us 
stress that in this preliminary study we have not included a transition region 
between the cold prominence and the hot corona (PCTR). The presence of a 
hotter plasma in the PCTR may somehow decrease the sensitivity of the Hell 
resonance lines to the Doppler effect as collisions will become more important 
in the formation processes of these lines. This will be investigated in a future 
work. 

The right panel of Fig. [1] indicates that there is no strong Doppler effect 
on the Call resonance lines, while we observe some Doppler brightening of the 
8542 A line (and indeed of the other two infrared lines at 8498 and 8662 A, not 
shown) at low temperature. 



3.2. Line profiles 

If spectroscopic observations of erupting prominences are available, then a com- 
parison between computed and observed line profiles can be made. We show in 
Figs. [2] and [3] the line profiles for the same helium and calcium lines considered 
in Fig. [T] at two different temperatures (solid line: 8000 K, dashed line: 15000 
K), at four different velocities (from top to bottom: 0, 80, 200, and 400kms~^). 

The Doppler dimming effect is well observed in the helium resonance lines 
at 584 A and 304 A as the radial velocity is increased (Fig. [2]) . We can observe 
asymmetries in the line profiles of these lines when the prominence plasma is 
moving radially, with some intensity enhancement which is especially visible in 
the red wing of the He I 584 A line at low temperature. This is explained as 
follows. The radiation emitted by the disk center in our code is represented by 
symmetrical line profiles. When the prominence is at rest, the incident radiation 
illuminating the prominence slab is symmetrical; however when the prominence 
plasma is moving radially the incident line profile becomes distorted and shifted 
towards the red. As we used the PRD standard approximation in our calcula- 
tions of the resonance lines of helium, the resulting line profiles are asymmetrical. 
This would not have been the case if we had considered complete redistribution 
in frequency (CRD) instead of PRD. The line asymmetry is more visible for 
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Figure 2. Line profiles for T = 8000 K (solid line) and T = 15 000 K (dashed 
line), with p = 0.1 dyncm"^, and L = 2000 km, at different velocities: 0, 80, 
200, and 400 km s^^ from top to bottom. Abscissa is AA in A and vertical axis 
is specific intensity in ergs~^ cm^^ sr^^ A~^. From left to right: Hel 584 A, 
Hell 304 A, and Hel 10830 A. 



the Hel 584 A line as its wings are fairly bright, while the wing intensities of 
the Hell 304 A line are low. Despite the fact that the line asymmetry increases 
with speed for both lines, it is more visible at low speeds (when the intensity in 
the wing is high enough). Finally, it is more pronounced at low temperatures 
because the scattering of the incident radiation is relatively more important as 
compared to collisional processes than it is at higher temperatures. 

Figure [3] shows that the intensities of the Ca H lin es are lower at 15 000 K 
than at 8000 K, an effect of CaH to CaHI ionization ([Gouttebroze &: Heinzel 
20021 ). 



4. Diagnostics of Thermodynamic Parameters 

We show here how the non-LTE radiative transfer calculations can help us to 
infer the thermodynamic properties of a prominence observed by the SUMER 
spectrometer on SOHO. This prominence was actually a rather quiet prominence 
and we have not included any velocity fields in these calculations. It was observed 
during the 13th MEDOC campaign held at IAS in June 2004. We select a few 
pixels in the SUMER slit which cut across the prominence and average the 
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Figure 3. Same as in Fig. [2] for Call K and Call 8542 A lines. 



line profiles there. We consider the line profiles of two H resonance hnes (Ly/3 
and Lye) and the He I resonance hne at 584 A. For the comparison between 
computed and observed line profiles we now include the presence of a PCTR. 
The tempera ture variation between the cold prominence core and the corona 
suggested by lAnzer fc Heinzell (|l999l ) has been adopted for this study. By a 
trial and error process we selected the temperature profile shown in Fig. SI 
The other prominence parameters are p = 0.023 dyncm"^, L = 194 km (total 
column mass 2.4 10~''gcm~^), H = 113722 km, and the microturbulent velocity 
^ = 18 kms~^ 

We obtain a very good agreement between the computed profiles (convo- 
luted with the SUMER instrumental profile) and the observed profiles, as shown 
in Fig. [5j It is worth noting that fitting hydrogen and helium resonance lines 
simultaneously places strong constraints on the parameter space, and it was not 
possible to find another set of parameters for the prominence that would be sig- 
nificantly different than what is given above and that would lead to a satisfactory 
fit of the observed profiles. 



5. Conclusions 

The non-LTE radiative transfer modeling that we are developing is a key tool 
for interpreting observations and constructing realistic prominence models. The 
combination of lines from hydrogen, helium, and calcium, places strong con- 
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Figure 4. Temperature as a function of the column mass in one half of the 
prominence slab. 




Figure 5. Observed (histograms) and computed (curves) lines profiles for 
Ly/3 (left), Lye (middle), and He I 584 A (right). 

straints on the models. Imaging and spectroscopy must be used for comparisons 
with calculations to determine thermodynamic parameters and velocities. The 
radial velocity determined from the comparison between observed and computed 
line profiles, in combination with line-of-sight velocities, should allow us to in- 
fer the full velocity vector of the prominence plasma. In a future work we 
will compare our model results with simultaneous observations of, e.g., Ha and 
Hell 304 A. 
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